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Abstract. We extend the ModeCode software of Mortonson, Peiris and Easther [1] to 
enable numerical computation of perturbations in K-inflation models, where the scalar 
field no longer has a canonical kinetic term. Focussing on models where the kinetic 
and potential terms can be separated into a sum, we compute slow-roll predictions for 
various models and use these to verify the numerical code. A Markov chain Monte 
Carlo analysis is then used to impose constraints from WMAP7 data on the addition 
of a term quadratic in the kinetic energy to the Lagrangian of simple chaotic inflation 
models. For a quadratic potential, the data do not discriminate against addition of such 
a term, while for a quartic (A0^) potential inclusion of such a term is actually favoured. 
Overall, constraints on such a term from present data are found to be extremely weak. 
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1 Introduction 

Observations, especially including those from the Wilkinson Microwave Anisotropy 
Probe (WMAP) [2-5], are beginning to impose useful constraints on inflationary cos- 
mologies. In particular, a number of papers have made detailed evaluations of con- 
straints on the simplest inflation models, featuring a single canonically-normalized 
scalar field with unknown potential V{(f)), see for example Refs. [4, 6-17]. 

Staying with the single-field paradigm, a more general scenario is available through 
the K-infiation paradigm. This retains minimal coupling of the scalar field to gravity, 
but allows the action to have an arbitrary dependence on the field's kinetic energy as 
well as on its value. This introduces new features, including a sound speed less than 
the speed of light which may enhance the non-gaussianity in the models. Our aim 
in this paper is to impose observational constraints on versions of this more general 
single-field scenario. 
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Our strategy is to modify the ModeCode program of Mortonson et al. [1], which 
solves the inflationary perturbation mode equations numerically and then interfaces to 
the CAMB [18] and CosmoMC [19] packages in order to compute the corresponding 
microwave anisotropies and compare to observational data. This entails a number of 
modifications to the way that ModeCode handles both the background (homogeneous) 
evolution equations and the perturbation equations. In this paper we will focus on 
the simplest case where the kinetic and potential terms remain sum-separable, and 
consider only simple forms for each. For the potential we will consider the simplest 
chaotic inflation models [20, 21], based on quadratic and quartic potentials. For the 
kinetic term, we will consider simple monomial and polynomial forms, in particular 
investigating constraints on addition of a term quadratic in the kinetic energy to the 
normal canonical form. Future work will explore more complicated forms. Modification 
of ModeCode to consider specifically the action corresponding to DBI inflation (see 
Refs. [22, 23]) has already been carried out in Ref. [24]. Comparison of K-inflation 
models to five-year WMAP data using slow-roll methods has been made in Ref. [25] . 



2 The K-inflation model 

The K-inflation model [26, 27] features a single scalar field with the action 

S= [ V^p{^,X)d'x, (2.1) 



where is the field value and X = (1/2)9^0(9^0. The function p will play the role 
of the pressure. A canonical scalar field has p{(f),X) = X — V{(j)) where V{(f)) is the 
potential. In this paper we will focus on models where the kinetic and potential terms 
can be written as a sum: 

p(0,X) = ir(X)-r(0), (2.2) 

where K{X) and V{(j)) are both arbitrary functions to be determined from data. 

Given the Lagrangian p{X, 0) of the considered model, we can obtain observable 
consequences by the following approach, closely following Garriga and Mukhanov [27]. 

2.1 Background field equations and sound speed 

We assume the usual Einstein equations and a spatially-flat Robertson-Walker metric 
with scale factor a and Hubble parameter H = a /a. We use the reduced Planck 
mass, defined by Mpj = (SttG)^^, to denote the strength of gravity throughout (with 
c = h = 1 as usual). 

First, without needing to consider gravity, we can obtain a relationship between 
the density p of the universe and its pressure p, 

p = ~3H{p + p), (2.3) 

which is known as the continuity equation, in which the energy-momentum tensor is 
characterised by the pressure p{X, 0) and the density 

p(X,0)=2X^-p (2.4) 
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of this universe. 

Second, now invoking a theory of gravity in the form of general relativity, we have 
the Friedmann equation for flat cosmologies 

Taking Eq. (2.3) together with Eq. (2.5), we can therefore study the evolution of the 
scale factor a and the field variable 0. 

One important quantity, called the 'sound speed', describes the properties of the 
(f) field. Regarding the field as a fiuid system, we can introduce as 

2 _ P,X _ P,X /r, 

where the comma denotes the partial derivative with respect to X. 
2.2 Perturbation mode equations 

The background state for the field is given by Eqs. (2.3) and (2.5). However, con- 
frontation of infiationary models with data requires us to evaluate the perturbations 
they predict. We will just quote the two significant equations, first derived by Garriga 
and Mukhanov [27]. The scalar and tensor perturbations are described by quantities 
V and u whose Fourier components in the longitudinal gauge satisfy 

^+(c,V ^)... = (2,7) 

^ ^ U _ <P^\ ^ „ (2.8) 

respectively. Here, r is the conformal time, subscript k denotes the momentum space, 
and the variable k relates to the comoving scale by A = 27r//c. The curvature pertur- 
bation C is related to v hy C, = v / z, while in a fiat universe the background variable z 
can be expressed as 

2 a^(P + P) 



In momentum space, the canonical quantisation modes Vk and Uk give two differ- 
ent classes of perturbations, the scalar perturbations measured by Vk and the tensor 
(gravitational wave) perturbations by Uk- From Eq. (2.7), we see that Cs/H plays the 
role of the 'sound horizon', which the k mode leaves by satisfying c^k = aH. The 
power spectrum of the scalar mode can be expressed by 

Vs{k) = , (2.10) 

ZTT^ Z 

while in Eq. (2.8) the tensor mode leaves the usual horizon by satisfying k = aH, with 
spectrum 

P.W^I^I^f. (2.11) 
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2.3 Power spectra and observables 

Solving Eq. (2.7) in the slow-roll approximation gives the expression [27]: 

1 



(2.12) 

CBk=aH 



The tensor mode Eq. (2.8) has its usual power spectrum 

2 



(2.13) 

k=aH 



'PI 

Having given the mode equations and their spectra, we can study the two most 
interesting variables, the spectral index Us and the tensor-scalar ratio r. These two 
quantities can be constrained by observations, such as the existing WMAP data or the 
forthcoming Planck satellite results. Their definitions are 

c? In P„ , , 

n,-l = — ~ -2e - 77 - s 2.14 
ct In A; 

nt = -rn^ ^ -2e 2.15 
a Ink 

r = § = 16ec,, (2.16) 

where the parameters e ,fj ,6 ,s are all small and defined as 

dlnH _ dine ^ din 6 dlnCs / «x 



dN ' ' dN ' dN ' dN 

through dN = Hdt, and we additionally give the tensor spectral index rit. Higher-order 
versions of these expressions, which we do not use here, have been obtained using the 
uniform approximation [28]. 

All these parameters are calculated at the time when mode k leaves its individual 
horizon. For the scalar mode k takes the value of aH/cs, while for tensor mode k takes 
aH. As a result, the relation of d/dN to d/dlnk for the scalars is 

dink , , 

— = 1-.-,. (2.18) 

Therefore Eq. (2.14) and Eq. (2.15) are then actually divided by Eq. (2.18), whereas 
s must be taken as zero when adjusting Eq. (2.15) as the sound speed does not enter 
its horizon-crossing expression. 

Note that fj in Eq. (2.17) is implicitly a function of the usual slow-roll parameters 
e = (M|[/2)(V^'/1/)^ and t] = Mp^V'/V and other smaller parameters, such as S.^ Here 
and throughout primes are derivatives with respect to 0. We will discuss its particular 
expression in later sections when investigating specific models. 



^The full expression is 

?7 = 2e - 2 ( 1 ■ ^^-^^ \ ^ . P- 



P.x J P.x ' 

given pressure p — p{X,(j)), where ,N indicates the derivative with respect to N. In our current 
consideration, the pressure has separable X and 4>, so the last term vanishes. The factor of 6 in the 
second term can be treated as —(1 + 6), where 6 = l/Cg. Therefore ry = 2e — (1 + 9)6. 
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3 Slow-roll predictions 



In this section, we use the slow-roll approximation to compute the spectral index ris 
and tensor-to-scalar ratio r for various models. These results are of interest in their 
own right as they give an indication of the properties of models that will be able to fit 
the data. Additionally, we will be able to use them to verify that our modifications to 
ModeCode have been implemented successfully. 

3.1 General prediction without specifying a potential 

We study models where the Lagrangian takes the form 

p(0,X) = ir„+iX"-l^(0), (3.1) 

where Kn+i are constants and n takes integer values.^ Under this kind of action, we 
can derive the field equation for the scalar field from Eq. (2.3) as 

^P,x + = -QnKn+iHX^ , (3.2) 

where the subscript , is the derivative with respect to the field. 

We now discuss what the models predict under the slow-roll scheme. Therefore 
by the usual consideration, we assume the density of the universe is dominated by the 
scalar potential, oc p ~ and replace all appearances of the Hubble parameter 
by V. In the field equation Eq. (3.2), we take the first (acceleration) term on the 
left-hand side to be much less than the second term. Hence it simplifies to 

y> ~ -QnKn+iHX'' . (3.3) 

Within the slow-roll assumption, we can obtain descriptions of the observable quantities 
by means of some parameters. From Eq. (2.17) we can identify e as under the slow- 
roll assumption, where 

IV' d) , , 

= (3.4) 

With <^ under the slow-roll assumption, for later discussion we give the second 
parameter rj^ which relates to the first and second-order derivatives of the potential. 

By the approximation above, we can find a simple relation between these parameters, 

ri VV" 

and also we can find the relation 

V = ^^v - ^ - Vv ■ (3-7) 



^With only a single kinetic term, the eoefficient Kn+i could be removed by rescaling </>, adjusting 
the potential, but for later comparison with cases with more than one term we keep it explicit. 
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From Eq. (3.3), for later evaluation we can write the ratio (p/H as 

| = -«(-)(^) (3-8) 

where 



nn-l \ l/(2n-l) 

a{n) = { ^^M2[^ . (3.9) 



riK, 



n+l 



Then 



1 



Y/2n \ l/(2n-l) / p^//(2n-l) \ l/(2n-l) 



ev = l^fTs^^J ; Vv=o^in) (^ ^,^,,,_, J • (3.10) 

Now we would like to apply the slow-roll approximation to evaluate the most 
important physical observables. The first is the scalar power spectrum Vs, 

Y5n-2 \ l/(2n-l) 



In this case, we can write the spectral index via Eq. (3.11) as 

ns-l = [^nr,^ - 2(5n - 2)e J . (3.12) 

To verify that Eqs. (3.11) and (3.12) are correct for standard inflation, we just need to 
set n = 1, and then we have the power spectrum 

V-^—^ (3 13) 

and its corresponding spectral index is ris — 1 = 2//, — 6e,. As a prediction of Eq. (3.12), 
for the simplest non-canonical inflation, in which the kinetic term has the form X^, we 
can obtain the power spectrum 



and its spectral index Ug — I = {4:1]^ — 16e,)/3. 

3.2 Predictions for specific models through e-folding 

Thus far we have obtained the formulae for the scalar power spectrum and its spectral 
index without specifying a particular potential type. In this subsection, we continue 
the discussion within a class of potential V{(f)) = Acf)^ in the Lagrangian Eq. (3.1), 
where A denotes the normalisation parameter, and m takes integer values.'^ 



■^Following footnote 2, for this form of potential a rescaling of to eliminate i^„+i in the single-term 
kinetic case would simply rcnormalize A, which anyway is to be fixed by the density perturbation am- 
plitude. Hence there is a perfect degeneracy between Kn+i and A and data can only fix a combination 
of them. 
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The e-folding number N, which measures how much inflation took place, is 



N = -jHdt. (3.15) 

After observable scales cross the horizon we have ~ 50 until inflation ends. First 
we evaluate the time variation H dt by the time-shift carried by the scalar field and 
finally transferred from the gradient change in the potential V{(f)) 

Hdt = ^d^=^—dV . (3.16) 

So from the definition of A^, and taking the above relation along with Eq. (3.8), we 
can derive the corresponding relation for A^ in terms of the scalar potential V{(j)) and 
its derivatives, rather than with itself, 

N^lj^^V. (3.1T) 

Under the assumption of slow-roll, by means of the slow-varying parameters, say 
Eqs. (3.4) and (3.5), we can write Eq. (3.17) as 

Here e (neglecting the subscript) is a function of the potential and its derivatives 
as well, which have already been denoted by Eqs. (3.4) and (3.5) in the slow-roll 
approximation. Connecting A^ with these two slowly-varying parameters leads to the 
expression for spectral index Ug. 

By taking potential V{(f)) = A(f)"^, we obtain 

p(n, m) 7(n, m) 

where 

^/ N m(n — l) + 2n , , , /^x 2?i/(2n-i) 

2n — I 

Then, connecting Eqs. (3.4), (3.5) and Eq. (3.19), we can get the relation between the 
slow-roll parameters and A^ 

ml m — 1 1 oT\ 

The final formula for the spectral index in Eq. (3.12) is therefore 

n,-l = -X(n,m)^, (3.22) 
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n = 1 


n = 2 


m = 2 


m = 4 


m = 2 


m = 4 


a{n) 






(3M|i/;r3)^/' 


(3M^,/ir3)^/' 


f3{n, m) 


2 


2 


2 


8/3 


j^n, m) 


Aa{l)A 


16a(l)Ai/2 


24/3«(2)A2/3 


28/3«(2)AV3 


X(n, m) 


2 


3 


2 


5/2 




7(i,2)iV^ 






y37?/y2iV/3)V2 


Us 


N^50 


0.96 


0.94 


0.96 


0.95 


r 


]V=50 


0.16 


0.32 


0.092 


0.14 



Table 1. Functions and observables for standard inflation and the simplest non-canonical 
inflation (NCI) model K{X) = K^X"^. Vs is in units of (Un'^M^^) and N is taken as 50. 



where 



X(n, m) 



■m{3n — 2) + 2n 



(3.23) 



■m[n — 1) + 2n 

We can separate the two dependencies in Eq. (3.23) so as to examine how the 
potential and kinetic energy terms contribute: 



(2n — l)m 
m(n — 1) + 2n 



1 + 



m 



(n , m > 1) 



(3.24) 



which indicates X > 1. We can see, from the above expression, that the spectral index 
1 has a simple relation with the power of the potential in A^, Eq. (3.19). In 



n. 



particular, for the quadratic potential m = 2, I is independent of n. Therefore, the 
inflation model driven by this potential will give a spectral index Ug = 0.96 regardless 
of the power of X in the kinetic term. Further, in terms of X{n,m) in Eq. (3.24), the 
scalar power spectrum Eq. (2.12) can be written as 



127r2M^j Cs 



7 



x-i 



N 



X- 1 



(3.25) 



The class of Lagrangian, Eq. (3.1), has sound speed Cg = l/\/2n — 1 which is 
independent of m. As for Eq. (3.25), given a known e-folding number N after inflation, 
the scalar power spectrum can be determined directly using the parameter set {n,m), 
which is the input argument to the two functions '-f{n,m) and I{n,m), while for the 
scalar spectral index we can just use one function X{n,m). Example values for the 
functions a(n), (3{n,m) and 7(77., m) with parameter set {n,m) can be seen in Table 1. 

From Eqs. (2.16), (3.21) and (3.22), we can obtain an expression relating the 
tensor-to-scalar ratio r and spectral index 



V2 



n 



1 m 



m(3n — 2) + 2n 



:i-ns 



(3.26) 



We list the observables in Table 1, and example values for and r can be seen in 
Figure 1. 
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Figure 1. Slow-roll predictions for the tensor-to-scalar ratio r and the scalar spectral index 
Us in standard inflation and the simplest NCI model where n = 2 in Eq. (3.1). 



For canonical inflation, the known relation rig — 1 
setting n = 1. For n = 2, we find 



4(m + 1) 1 
m + 4 



(m + 2)/2N is recovered by 



(3.27) 



Figure 2 shows the spectral index as a function of m for several n values; when n > 1 
the spectral index asymptotes to 1 — in the limit of large potential power-law m, 
unlike the canonical case where 1 — ris grows linearly with m and can be large. 



4 ModeCode for K-inflation 

For single-field canonically-normalized inflation models there are many numerical tools 
which can calculate the primordial power spectra, as well as other characteristics such 
as the bispectrum and trispectrum of non-gaussianities. Also many of these tools have 
been interfaced with MCMC codes such as CosmoMC so as to explore the likelihood 
and carry out parameter estimation. ModeCode [1] is an example of such a programme, 
which had also recently been interfaced to the MultiNest model selection code [29, 30]. 
Other codes to numerically solve the inflationary mode equations have been described 
in Refs. [31-36]. 

To study non-canonical models, modifications must of course be made to the code, 
generalizing both the background and perturbation equations. An example already in 
the literature is an enhancement to consider the Lagrangian describing DBI inflation, 
made in Ref. [24], which is motivated from string scenarios [22, 23] and has been 
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m (Power of Potential) 

Figure 2. Slow-roll predictions for the spectral index with different kinetic power-law n and 
potential power-law m. This plot assumes the pivot corresponds to = 50. 

widely studied in the literature. In this paper we consider a different extension to non- 
canonical inflation models, at this stage restricted to models where the kinetic energy 
K{X) and the potential V{(f)) are sum-separable. 

4.1 Brief description of ModeCode 

ModeCode is characterised by free parameters describing the inflationary potential for 
canonical single-field inflation, and solves the inflationary mode equations numerically, 
bypassing the slow-roll approximation. It computes the cosmic microwave background 
angular power spectra and performs a likelihood analysis and parameter estimation by 
interfacing with CAMB [18] and CosmoMC [19]. 

4.2 Modifications needed for K-infiation 

We implement ModeCode under our umbrella Kinetic Module Companions (KMC) . 

To perform our analysis for NCI models, we build up a full system for initialising the 
background equations, avoiding slow-roll or any other kind of approximation beyond 
linear perturbation theory. The parameters which can be explored and the methodol- 
ogy in KMC are as follows. 

• Parameters describing the form of the kinetic energy, for example a Taylor ex- 
pansion of K{X) about X = which we have currently implemented up to sixth 
order (though in the present paper we will only consider up to quadratic order). 

• Parameters describing the inflationary potential, for instance the potential can 
take a polynomial form or be a Taylor series. 
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slow-roll predictions and numerical results 
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Figure 3. Slow-roll predictions and numerical results for standard canonical inflation and the 
simplest NCI model [K^X"^ — Acj/^) , showing they are consistent. The matter power spectrum 

amplitude constrains the combination AK^^^ (from the scaling argument of footnote 3) and 
the predictions for and r are independent of this. The value of K^, does determine the 
e-folding value corresponding to observable scales, but this is not fixed by observations. 



We use eigenvalue methods to get the real solutions needed by the background equa- 
tions, Eqs. (2.4) and (3.3), as well as perturbed equations. Equations to solve simul- 
taneously for and can be explicitly obtained. 

4.3 Comparison tests 

4.3.1 Recovery of ModeCode results for K{X) = X 

Before running the extended functions, we check we can recover the results of Mode- 
Code with our KMC system. The outputs that ModeCode generates are recovered 
either by setting our flag use_kinetic=T and setting K{X) = X, so that KMC will 
perform its intrinsic functions for models that can be executed by ModeCode, or just 
by switching off the KMC functional system by flag use_kinetic=F, leaving KMC to 
function as the normal ModeCode. We have confirmed that the ModeCode results are 
then precisely recovered under either method. 

4.3.2 Recovery of slow-roll results 

Now we compare the slow-roll predictions of the previous section with the numerical 
results for two types of inflation model. One is the standard inflationary model, with 
canonical kinetic energy and Lagrangian p{(f), X) = X — V{(j)). The other is the simplest 
non-canonical inflation model, where the kinetic energy takes the form X"^ in the 
Lagrangian p{(f),X) = K^X"^ — V{(j)). Figure 3 shows that the recovery of the slow- 
roll results is very accurate. It is not expected to be absolutely precise because the 
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slow-roll approximation is not perfect, for instance leading to an offset in identification 
of the = 50 point as well as neglecting higher-order corrections to perturbation 
observables. 

5 Parameter explorations with MCMC methodology 

5.1 Global settings and initial conditions 

We now proceed to explore the parameter space of a particular non-canonical type of 
inflationary model, by means of Markov Chain Monte Carlo (MCMC) methodology. 
We use the WMAP 7-year data version4 ('WMAP7"), and our simulations take 12 
chains for each model. We set the pivot scale to /cpivot = 0.05 Mpc^^. We aim to select 
an initial field value 0init which corresponds to 70 e-foldings from the end of inflation, 
estimated analytically by assuming a single power-law term dominates K{X); if this 
approximation proves inaccurate it gets adjusted by the numerical code. We then must 
choose a consistent initial field velocity cf) j^mit, as mentioned in Section 4.2, a task which 
has been solved by eigenvalue methods in our KMC numerical modules. 

Having developed the KMC code, our aim is to investigate what type of NCI 
models are supported by observational data. Using the MCMC method, we will per- 
form a likelihood analysis and parameter exploration for some particular non-canonical 
inflation models. 

5.2 Choice of models 

In this article, for our numerical work we focus on a particular choice of kinetic term 
which adds a quadratic term in X to the usual linear one. Investigation of more 
complex models will be made in future work. Hence our considered NCI model is 

0) = K2X + K^X^ — V{(j)) with positive,'^ and we will additionally assume 
the potential to have a single polynomial term V{(j)) = Acjf^ with m = 2 or 4, giving 
a large- field model. The field (f) can always be rescaled to set K2 = 1., and the new 
term with coefficient can be considered as the first correction term in a Taylor 
expansion of a general K[X) that reduces to a canonical form in the limit X — )■ 0. 
While such a model is not particularly realistic, it has the benefit of simplicity and it 
is interesting to ask whether present data can say anything about the possible values 
of such a correction. Since in slow-roll inflation X will be small, we can immediately 
anticipate that any constraint on will be very weak, allowing values much greater 
than one before this correction term could significantly modify the canonical term. 

Our free parameters are therefore the power of the potential, which we fix for 
each investigation, and the values of the amplitude of the potential and the coefficient 
K^. Additionally, the value of A^pivot corresponding to A;pivot can vary as in the orig- 
inal ModeCode. The well-measured amplitude of perturbations will accurately fix a 
combination of these parameters. 

^Negative appears possible in principle, provided X is not too large, but gives models that 
can have phantom behaviour [w < —1) from an overall negative kinetic term, which can be expected 
to cause instabilities. We note also that the quadratic approximation to the DBI model features a 
positive coefficient. 



- 12 - 



Models with fixed log-n 



-10 



Models with fixed log-io X = 





STD (X-mV/2) 




NCI (log k3=13) 




NCI (log k3-12) 




NCI (log k3=11) 




^y"^ NCI (log k3=10) 




\ NCI (log k3=9) 




N^''S. NCI (log k3-8) 







200 400 600 800 1000 1200 1400 1600 1800 2000 

(a) plot for quadratic potential 
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(b) plot for quartic chaotic potential 



Figure 4. The effect of different values of on the final power spectrum. All numbers in 
the legend are base- 10 logarithms, with the spectrum in units of ji¥i?. 
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Figure 5. The power spectra for various combinations of parameters. For the NCI models, 
the first number in the key is the exponent of m? (left panel) or A (right panel), while the 
second number is the exponent of K^. 



5.3 Overview of effects from additional kinetic terms 

Before we examine the results from MCMC, we look at the influence of the extra term, 
i^aX^ in the kinetic function, on the final power spectrum. The shape of the power 
spectrum is controlled by a combination of A and in our considered NCI models. 
In the left panel of Figure 4, we take logm^ = —10 for potential m^0^/2; introducing 
has negligible effect for < 10^°, above which the spectrum starts to decrease 
as the quadratic kinetic term becomes dominant. In the right panel we see a similar 
result for with A = 10~^^. Figure 5 shows the spectra for some parameter values 
chosen so that the spectral amplitude is close to the observed value. 



5.4 Interpretation of MCMC explorations 

After finishing 12 chains for each model, we have obtained a Gelman-Rubin conver- 
gence statistic R — 1 = 0.16 and 0.012 for the eigenvalues of the covariance matrix in 
models with quadratic and quartic potentials respectively. The prior ranges and max- 
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Figure 6. Parameter constraints for NCI model with quadratic potential. Left: Constraints 
on logm^ and against A'^pivot for WMAP7 data. The countours are 68% (inner) and 
95% (outer) confidence levels, while the colour scale shows the sample mean likelihood in 
bins. Right: one-dimensional posterior distributions for the parameters (solid) and the mean 
likelihoods (dashed). 



imum likelihood values are given in Table 2, and the posterior distributions in Figures 
6 and 7. 

The parameter of principal interest in each case is K^. In the quadratic potential 
case, this parameter turns out to be completely unconstrained by the data. This is to 
be expected, as the quadratic potential gives acceptable observables when the kinetic 
term is either X or X^, which are the limits of small and large K^. The MCMC results 
show that the fit to data remains acceptable right to the largest values of that we 
permit. We clearly see the two limiting behaviours of domination by either the X or 
X"^ term; for example in the 2D m^-K^ constraint plot the former region has constant 
logm^ ~ ~10, while the latter has i^s oc as implied by taking constant Vs in the 
third column of Table 1 to obtain the observed amplitude. The bimodal likelihood of 
r is caused by the different values of this parameter in the two regimes. We also see 
a very mildly enhanced likelihood in the transition regime ~ 10^^, but all values 
of are acceptable Perhaps surprisingly, then, present data can say nothing about 



Models 


Priors 






-2 In £ml 


(NCI) 


log A log 








(1,2; 2) 


(-16,-4) (0,20) 


0.965 


0.080 


7469.8 


(1,2; 4) 


(-18,-4) (0,20) 


0.957 


0.115 


7471.8 



Table 2. Priors for the parameters, log A and logA'3, and the maximum likelihood (ML) 
values for ris and r from WMAP7. NCI models with single-term potentials provide a red tilt, 
ns < 1, and a detectable tensor-to-scalar ratio r ~ 0.1. 
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Figure 7. Parameter constraints for NCI model with single term quartic potential. This 
figure uses the same convention as in figure 6. 



the ampUtude of a quadratic kinetic term added to the normal canonical one for this 
potential, and as this is a potential known to fit the data well in the canonical case we 
can conclude that more generally a quadratic correction term cannot be constrained 
directly from data. 

The quartic potential case shown in Figure 7 is more interesting. In the standard 
cosmology, — >■ 0, the quartic potential is quite disfavoured by WMAP7. As we 
would anticipate from the slow-roll results of Section 3, the situation for this potential 
actually improves in the limit of large K3, as both ng and r move towards the scale- 
invariant values. Once more the MCMC results show that the fit remains good as the 
quadratic kinetic term becomes dominant and we find no observational upper limit on 
for this potential. This time the distribution for r is unimodal, as only the X"^ 
domination regime contributes to the posterior. There is a plausible non-zero lower 
limit on K^, though the numerical value of such a limit will be quite prior-dependent. 
Hence, incorporation of an X"^ term is a method of salvaging the quartic model, though 
the large value of K^, in Planck units, that is required is unattractive. 

6 Conclusions 

In this paper we introduced a numerical solver Kinetic Module Companions (KMC), 

an extension to ModeCode for a class of non-canonical infiation (NCI) models. In this 
article we have used our code to investigate some simple non-canonical models, which 
have up to two terms in the kinetic energy and a monomial potential, in order to test 
the validity of the code and provide some initial scientific results. We found that these 
models are well able to fit current data, including in the quartic potential case pro- 
vided the quadratic kinetic term dominates. This is compatible with slow-roll results 
we obtained in the case of a single kinetic term of power-law form. 



- 15 - 



As a specific application of tlie code, we studied tlie introduction of a quadratic 
correction K^X^ to tlie normal canonical kinetic term X, with the goal of constraining 
the coefficient K^^. In practice, however, turns out to be unconstrained by data, 
and indeed the inclusion of a large quadratic term can even improve the fit to WMAP7 
data, for instance for a quartic potential which, in the canonical case, is under severe 
pressure from observations. Accordingly present data allow no leverage whatsoever 
on radical deviations from the canonical case. In future work we plan a much more 
comprehensive investigation of possible kinetic and potential forms. 

A longer-term objective in this area may be to extend ModeCode to yet more 
complex forms of single scalar-field action, such as the Galileon [38] or indeed the 
Horndeski action [39] which is the most general scalar-tensor theory yielding second- 
order equations of motion. However the many functional degrees of freedom of such 
actions will no doubt lead to considerable degeneracies given the relatively limited 
amount of observational information available, which essentially amounts to only a 
couple of numbers at present. Hence, as we have found here even for the simplest 
separable K-infiation case, one is likely to need considerable guidance from theory as 
well as from observations in assessing whether the most general paradigms are useful. 

Noted added: Shortly after we put this article on the arXiv, an independent paper 
[37] was submitted which contains results that partially overlap with ours. 
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